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Abstract. We present a simple and fast algorithm for the computation of the 
coefficients of the expansion of a function /(cos u) in ultraspherical (Gegen- 
bauer) polynomials. We prove that these coefficients coincide with the Fourier 
coefficients of an Abel-type transform of the function /(cos u). This allows 
us to fully exploit the computational efficiency of the Fast Fourier Transform, 
computing the first N ultraspherical coefficients in just 0(A r log 2 A r ) opera- 
tions. 



1. Introduction 

In this paper we investigate the generalization to ultraspherical polynomials (also 
known as Gegenbauer polynomials) of the results obtained in a previous paper [5] 
regarding the efficient computation of the coefficients of Legcndrc expansions. 

Ultraspherical expansions play a relevant role in various subjects of applied and 
computational mathematics. These expansions have been used successfully for the 
solution of linear |8] and nonlinear j6[ [7j 115j differential equations, of integral equa- 
tion [Jj)] , and in spectral methods for partial differential equations [21 [Jj5] . Gegen- 
bauer filtering (for the suppression of the Gibbs phenomenon) [TQl [Til HI] has been 
proved to provide an exponentially convergent approximation (in the maximum 
norm) of a piecewise analytic function, starting from the Fourier partial sum of the 
function itself. The related Gegenbauer reconstruction method has found natural 
application in the image segmentation problem, in particular of MRI images [H [Jj5] . 
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One of the limits of the application of ultraspherical expansions is the high cost 
of computing the expansion coefficients, this question becoming particularly critical 
in the case of multivariate functions. 

In this paper we give an efficient procedure for computing these coefficients. We 
first obtain a Dirichlet-Murphy-type integral representation of the ultraspherical 
polynomials pi^ (cos u) of degree n and order d. Then we prove that the coefficients 
of the ultraspherical expansion of a function /(cosu) (u € [0, 7r]) coincide with the 
Fourier coefficients (restricted to nonnegative index) of an Abel-type transform of 
the function /. 

These results produce straightforwardly an algorithm for the computation of 
ultraspherical coefficients which is very simple and very fast. In fact, the first N 
ultraspherical coefficients are obtained in only 0(7Vlog 2 N) operations by a single 
Fast Fourier Transform of the Abel-type integral function, the latter being easily 
computable by standard quadrature techniques [14]. 

Finally, in the Appendix the results obtained for the ultraspherical polynomials 
are given also in terms of Gegenbauer polynomials Cn (x) with x G [—1,1]. 



2. Connection between Ultraspherical expansions and Fourier series 

In dimension d, the expansion in ultraspherical polynomials of a function / = 
/(cosh) (u € [0,7r]) reads: 
(2-1) 

/(cosh) = Y F(n + f~ 2) (n + ^) 4 rf > P ^(cosu), 

2 (d-2) ^r^A. n \ \ 2 ; 

the coefficients a n being defined by [I]: 

(d-l) 

(2.2) = f /(cosh) P W(cosu) (sinn)( d - 2 ) du (n e N ), 

1 { — ) Jo 

where P„^(cosii) denotes the ultraspherical polynomial of degree n and order d, 
which is defined by the integral representation [SJ [TB] : 

v^r(^) 

Our goal now is to show that the coefficients affl are the Fourier coefficients of 
a suitable Abel-type transform of /. To this end, we first prove the following 
proposition (see also [4]). 

Proposition 1. The following integral representation of the ultraspherical polyno- 
mials P„ (cosm) holds: 
(2.4) 

pM>(oob«) = T{ i 2 \ 1 / e^+^ [2(cos M -cost)]^ dt. 

Proof. In the integral representation (|2.3|) of Pn d ^ (cos u) substitute to r\ the complex 
integration variable r defined by 

(2.5) e lr = cos u + isinu cost;. 



(2.3) P„ (rf) (cosu) = 7=^^ / (cosu+isinucosyy)" (sin??) (d - 3) dr? (n G N ). 
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It can be checked that 

(2.6) 2e ir (cosr - cosw) = (e iT - e in )(e ir - e~ iu ) = sin 2 u sin 2 ??. 

Now, since e lr dr = — sin u sin ij dij, the integrand on the r.h.s. of Eq. (|2.3p can be 
written as follows: 

- (sinu)-( d - 3 V (n+1)T \(e iT - e iu )(e iT - e"' 11 )] ^ dr 
= -(sinw)- (d - 3) e i( " + ^ )T [2(cosr- cosw)]~ dr. 

In order to determine the integration path, consider an intermediate step where e 1T 
is chosen as the integration variable; the original path (corresponding to n 6 [0,7r]) 
is the (oriented) linear segment Sq(u) starting at e™ and ending at e~ lu . Since (as 
shown by (|2.7p ) the integrand is an analytic function of e 1T in the disk |e lr | < 1 
(since n € N), the integration path 5q(u) can be replaced by the circular path 
5+(u) = {e 1T ; r = t, u ^ t ^ 2-7T — u} (see Fig. [TJ. Moreover, by using the fact that 
[(e iT - e iu )(e iT - e - *")]^ is positive for e ir e <5 (u) U M and therefore at e iT = e i7r , 
we conclude from the left equality in (12 .6[) that in the r.h.s. of (|2.7I) the following 
specification holds (for r = t; u ^ £ 2tt — u): 

(2.8) [2(cost - cosm)]^ = (-i) d ~ 4 [2(cosw- cost)] 5 * 4 . 

Finally, by taking into account the latter expression, the integral representation 
(|2.3|) can then be replaced by the integral representation (|2.4[) . □ 

We can now prove the following theorem (see also [4]). 

Theorem 2. The ultraspherical coefficients {a^ }*^L (see (|2.2p ) coincide with the 
Fourier coefficients (restricted to the nonnegative integer) of the form: 

(2.9) agO - fW(t)e int dt (n G No), 
(2.10) 



= [-ie(t)] d ~ 2 ZZ^ I /(00s u) [2(cosu-cosi)]^sinudu, 

2~ 



r(^) 



e(i) feeing i/ie sign function. 

Proof. By plugging representation (12.4[) into formula (|2.2[) we have: 
(2.11) 

alf> = (~i) d - 2 . / dii/(cosu)sinu / e i( " + ^)* [2(cosu - cosi)]^ di 

r(— 2— ) Jo Ju 

= (-i) d ~ 2 r(rf _ 2) |y dte i(n+ ^)*y dw /(cos u) sin u [2(cosw- cost)] t 4 

f 27r d 2 r 2 ™-t d _ 4 

+ / di e^ nH a - )*/ dw /(cosu) sinu [2(cosw — cosi)]~ 
Jtt Jo 

Now, changing the variables (t, u) —> (t + 2tt, u) in the second term inside the 
parentheses on the rightmost side of (|2.11[) , the latter becomes: 

pO _ p-t 

(2.12) e ™ {d - 2) / dte i( "+V)*/ du /(cos u) sin u [2(cosu - cost)] 
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Finally, from (gHH) and (12TT2]) we have: 
(2.13) 



4? = (-i) 



<Z-2 



2tt— 



df 



e 1 2 * / dw /(cosu) sinu [2(cosu — cost)] - 



,0 




/ dt e mt 




/ — 7T 


Jo 



du / (cos u) sin u [2 (cos u — cos t)\ 2 



e int dt, 

with J {d) {t) given by (|2~TUj) . 



□ 



It is easy to check that the 27r-periodic f^ d \t) function enjoys the following sym- 
metry properties: 

(2.14) fW(t) = (-i) d e 1 ( d - 2 ) t /W(-t) (ten), 

and, consequently, 

(2.15) aW=(-l) d a^ n+d -2) i^eZ). 
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Appendix A. From ultraspherical to Gegenbauer polynomials 

A.l. The Gegenbauer polynomials . For the convenience of the reader we 
summarize hereafter the main properties of the Gegenbauer polynomials Cn(x), 
which are strictly related to the ultraspherical polynomials P n (x) that we intro- 
duced in Section [2j The following formulae are listed in Ref. [20] and can also be 
found in refs. [Q Chapter 22] and [TTJ Q7] . 

The Gegenbauer polynomials of order A can be defined in terms of their generating 
function: 



(A.l) ^Ci A )(x)r = (l-2 a; t + t 2 ) 



— A 



From (|A.1|) we see that the Legendre polynomials P n (x) are the particular case of 
Gegenbauer polynomials with A = i, i.e.: Cn 2 \x) = P n (x). 

They satisfy the recurrence relation: 
(A.2) 



nC^(x) = 2(n + A - l)xC£\(x) - (n + 2A- 2)cj A) 2 (x), 



C (A) (a;) = l, c} X) (x) = 2Xx. 
In terms of Gaussian hypergeometric function they can be written as: 
(A.3) cWfc) = 2l ~ 2X ^Zt + 2A) ^ U 2A + n; A ' 1 ■ 1 ~ ^ 



n\ T(A) \ 2' 2 

They can be written explicitly as 

L"/2J TV\-L _ h\ 

(A.4) g. w W=E(-')' M (i. + ayrw ' to »-"- 

and have the following integral representation: 

(A.5) Cj A) (cosu) = , F V ^ N10 — -/ (cosu + isinucos77) n (sin?7) 2A_1 d7i. 

n! [r(A)] 2 J 

The Gegenbauer polynomials can be computed by the Rodrigues formula: 

(A 6) C W(x) = ( - 2) " r(n + A)r( " + 2A) (l - x*Y X+ -— \(l ~ x 2 ) n+X -i 
[ > " { > n! r(A)r(2n + 2A) 1 j dx™ l l Xj 

which follows by induction from 

(A.7) A C W(x)=2ACi A l 1) ( a ;). 

For fixed A, the Gegenbauer polynomials are orthogonal on the interval [—1, 1] with 
respect to the weight function: 

(A.8) w w (x) = {l-x 2 ) x -^, 

that is, for n^m: 

(A.9) £ C„W(x) (1 - x 2 ) x -i dx = 0, 
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and are normalized by: 

J -i L n\ (n + A) [r(A)] 

A. 2. Relation between Gegenbauer polynomials and ultraspherical 

polynomials Pjf 1 ' • Comparing the integral representations (|2.4|) and (|A.5|) . it is 

easily seen that the relation between ultraspherical P„ ^(cosw) and Gegenbauer 
polynomials ci A ^ (cos it) is: 

(A.lla) CWfooBt*) = P„ (2A+2) (cos W ), 

(A.llb) Pj d )(cosu) = r ^ r |V-2) 

Then, rephrasing Proposition [1] in terms of Gegenbauer polynomials, from (j2.4[) 
and (|A.11[) it follows that the polynomials Cn^ (cos u) (u € [0, 7r]) have the following 
integral representation: 
(A.12) 

°" ^ co& ^- ^n!r(2A)r(A) (sinu) 2 *" 1 6 ^ C ° &M C ° SfjJ dt " 

A. 3. Expansions in Gegenbauer polynomials . The Gegenbauer expan- 
sion of a function f(x), defined in x G [—1,1], reads [TTj : 

oo 

(A.13) f(x) = £ /n (A) Cn (A) W (* e [-i, i]; A > o), 

where the Gegenbauer coefficients are given by 
(A.14) 

f^= fSl^t^l t /(»)^(») d-^)(W)d, („ 6 No). 
y/TT T(n + 2A) T(A + 5 ) J-i 

Now, from (|A.14|) . (IA.11[) and using definition (|2.2[) of the coefficients al we see 
that the Gegenbauer coefficients arc related to the ultraspherical coefficients by 

We can now rewrite Theorem [2] in terms of Gegenbauer coefficients fn X ^ ■ 

Theorem 2'. The coefficients {(fffl}'£L Q , defined as (jffl == ^" +x ^ ( see (|A.14[) ) co- 
incide with the following Fourier coefficients (restricted to the nonnegative integer): 

r(\) p 7T 

(A.16) ^ = = / <f>^(t)e int dt (n e N ), 

(n + A) 

(A.17) (A) (t) = [ ~ 1£(f ^ 2A e iAt / f(x) [2 (x - cosi)]^ 1 da?, 

n J cost 

e(t) being the sign function. 
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